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Abstract. 

Using density-functional calculations, we study the effect of sp 3 -type defects created 
by different covalent functionalizations on the electronic and magnetic properties of 
graphene. We find that the induced magnetic properties are universal, in the sense 
that they are largely independent on the particular adsorbates considered. When a 
weakly-polar single covalent bond is established with the layer, a local spin-moment 
of 1.0 fiB always appears in graphene. This effect is similar to that of H adsorption, 
which saturates one p z orbital in the carbon layer. The magnetic couplings between 
the adsorbates show a strong dependence on the graphene sublattice of chemisorption. 
Molecules adsorbed at the same sublattice couple ferromagnetically, with an exchange 
interaction that decays very slowly with distance, while no magnetism is found 
for adsorbates at opposite sublattices. Similar magnetic properties are obtained if 
£S) . several p z orbitals are saturated simultaneously by the adsorption of a large molecule. 

C*~) " These results might open new routes to engineer the magnetic properties of graphene 

derivatives by chemical means. 



PACS numbers: 73.22.Pr, 73.20.Hb, 75.70.Ak, 75.75.-c 
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1. Introduction 

Graphene has exceptional electronic properties with the potential to give rise to 
interesting new applications [TJ [2J |3|. In particular, graphene is very attractive for 
spintronics [4]. This is mainly due to the very long spin relaxation and decoherence 
times in these materials [5j [6]. There have been also numerous claims of the 
existence of intrinsic magnetism in graphene associated with the presence of structural 
defects [74 El El UM CEU CE21 EES] , edge states pU [15], and partial hydrogenation of the 
layer. [161 El EEB] These observations suggest that graphene, besides being an ideal 
material for spintronic circuitry, could also be used in active spintronic devices. An 
interesting example is provided by graphene nanoribbons. The energies of the polarized 
states localized at the edges of zigzag nanoribbons or nanotubes can be controlled by 
external electric fields, and this effect can be used to design a spin-filtering device. [T^[[T9] 
However, in spite of the growing interest, to date there has not been a conclusive 
experimental confirmation of this intrisic magnetism in graphenic materials. While 
some experimental groups claim to have found clear indications of ferromagnetism 
in irradiated graphite [7J, graphene prepared from reduced graphene oxide [20] and 
hydrogenated graphene [21] . these results are challenged by other groups with alternative 
experimental evidences [221 123] . 

One of the main difficulties for the observation of ferromagnetism in defective 
and hydrogenated samples might be related to the bipartite lattice of graphene. It 
is theoretically well known, at least within a simplified description based on a ir- 
tight-binding Hubbard model, that defects created in different sublattices tend to 
couple antiferromagnetically [2U [101 [251 [26] . Since, for example, defects created by 
irradiation are expected to be randomly distributed, this would prevent the formation 
of ferromagnetic order. Therefore, to generate ferromagnetism in graphene and other 
graphenic nanostructures, it is desirable to create all the defects or to adsorb all the 
hydrogen atoms in one of the graphene sublattices. [18] Unfortunately, technically this 
might be very difficult to achieve. For this reason, it is interesting to look for alternative 
ways to induce magnetism in the carbon layer. In the present work we analyze, using 
first-principles electronic structure calculations, the role of different sp 3 -type defects 
created by covalent functionalization as a source of magnetism in graphene. Covalent 
functionalization has not been frequently investigated in this context, even though it 
has been used for a long time in the chemistry of carbon nanostructures and now in 
graphene [271 [281 [29]. We find that, under very general conditions, any molecule attached 
to the carbon layer through a weakly-polar single bond produces an effect similar to 
that of hydrogen adsorption, i.e., induces a significant spin moment in the system. 
In particular, we explicitly show that the induced magnetism is largely independent 
on the size of the adsorbed molecule and even on its chemical or biological activity. 
Regarding the coupling of the induced moments, we show that identical results are 
obtained by the simultaneous adsorption of several molecules or by the adsorption of a 
larger molecule that creates several bonds with the layer. This simple fact may open 
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new and interesting routes to tune the magnetism of graphene. In principle, molecules 
can be produced with the appropriate structural and chemical characteristics to create 
bonds with the layer according to predefined patterns. Thus, the powerful techniques 
of organic and surface chemistry could be applied to synthesize magnetic derivatives of 
graphene that behave according to well studied theoretical models [211 HOI [251 EE] - This 
is particularly attractive due to the recent sucessful synthesys of different graphene 
derivatives using surface chemical routes [301 EE]- Thus, the synthesis of carbon 
nanostructures with functional groups at predefined positions, starting from previously 
functionalized monomers, seems plausible nowadays. Our study deals with the magnetic 
properties of such nanostructures. 

In more detail, our calculations predict that, when a single C— C covalent bond 
is established between an adsorbate and graphene, a spin moment of 1.0 \i& is always 
induced in the system. The size and spatial distributions of these moment are nearly 
independent of the particular adsorbate. We explicitly show that this effect occurs for 
a wide class of organic and inorganic molecules with different biological and chemical 
activity, e.g. alkanes[32]. polymers [33| 134"! 135] , diazonium salts[36l E3 EH], aryl and 
alkyl radicals [391 HO] , nucleobasesjHJ 02], amide and amine groups [13]. glucose, [S] and 
organic acids [38l |4"5] . While a similar behavior is observed for hydrogen adsorption, we 
show for other adsorbates (NH 2 [16], OH [15] and F [Ul 08]) that the induced spin- 
moment decreases with increasing bond-polarity. Moments induced by adsorption in 
the same sublattice align ferromagnetically, with an exchange coupling that falls off 
very slowly with the distance between adsorbates (~ r~( 1+e \ e ~ 0.20). In contrast, for 
molecules adsorbed on opposite sublattices no magnetic solutions could be stabilized. 

2. Methods 

Our first-principles electronic structure calculations use density functional theory [19] 
as implemented in the SIESTA code [50]. We use the generalized gradient 
approximation [5]], norm-conserving pseudopotentials[52] and a basis set of numerical 
atomic orbitals [50]. A double-C polarized (DZP) basis set has been used for the 
calculation of the magnetic moments and electronic band structures of all our systems. 
We have checked that the relaxed structures with a DZP basis are almost identical 
to those obtained using a double-^ (DZ) basis. Therefore, we have used the smaller 
DZ basis for the relaxations of systems containing more than 100 atoms. The force 
tolerance for structural optimizations was 0.04 eV/A. The integration over the Brillouin 
zone was performed using a well converged /c-sampling equivalent to 136 x 136 x 1 k- 
points for the unit cell of graphene. The fineness of the real-space grid used by SIESTA 
was equivalent to a 150 Ry plane- wave cutoff. We have repeated some calculations 
with the VASP code [531 El] using a well-converged plane- wave cutoff energy of 400 eV 
combined with the pro jected-augmented- wave (PAW) method. Other computational 
details were as previously described. The use of PAW potentials allows to check the 
inherent limitations of the norm-conserving pseudopotentials. The two methods provide 




Figure 1. (Color online) (a) Schematic representation of the adsorption geometry 
of our systems and the nomenclature used in Table [1] X denotes here the atom 
directly bonded to graphene. Panels (b)-(e) show the isosurfaces (±0.019 e~/Bohr 3 ) 
of the magnetization density induced on graphene by functionalization using adenine 
derivatives [42 , CH 3 , PMMA and PTFE. Majority and minority spin densities 
correspond respectively to light and dark surfaces, which localize in different graphene 
sublattices and show a slow decay with distance in all the cases. Panels (f) and (g) 
show the spin polarized band structures of a 8x8 graphene supercell with a single 
adenine derivative and a CH3 molecule chemisorbed on top of one carbon atom, 
respectively. The black and red (gray) lines denote the majority and minority spin 
bands, respectively. The Fermi energy is set to zero. 

almost identical results. We have checked that, when a sufficiently dense k-sampling 
is used, our results are very weakly dependent on the supercell size. In particular, a 
spin moment of 1 //b is always associated with the formation of single C-C bonds with 
graphene, as well as with H adsorption. All the calculations presented in the text used 
a 8x8 supercell of graphene. 
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3. Results and discussion 

We first calculated the adsorption geometry of several molecules that can establish single 
bonds with the graphene layer. Most of them are attached to graphene through a single 
C— C covalent bond (denoted by X=C in Table [TJ X being the atom that bonds to 
graphene). We have also considered other adsorbates (H, NH2, OH and F) that form 
heteropolar bonds with graphene. In all cases the adsorbate sits directly on top of one 
of the C atoms of the layer [see the scheme in Fig. [JJ (a)]. This central C atom (Ci) 
moves vertically and approximately keeps a symmetric three-fold structure respect to its 
neighbors (C2). Table [TJ describes adsorption geometries of the different molecules. The 
angles Qx and 02 are in the range 102.1°— 107.0° and 111.8°— 115.7°, respectively. These 
values are half-way between those corresponding to a sp 2 hybridization of the central 
Ci atom (Qx = 90° and 62 = 120°), and those expected for a sp 3 local hybridization 
(Ox = ©2 = 109.5°). Thus, the observed geometries can be understood by a slight 
modification of the sp 2 hybridization of C\ that gains certain sp 3 character [55]. The 
largest sp 3 character among the studied molecules corresponds to the adsorption of 
nitrobenzene (C6H4NO2), anisole (C6H4OCH3), PMMA, and polystyrene groups. The 
smallest local distortion of the carbon layer corresponds to F adsorption, reflecting a 
larger ionic character of the interaction. 

The main results for the magnetism of all the studied adsorbates are presented 
in Fig. [1] and Table [TJ We first focus on the adsorbates attached to the layer by a 
single C— C bond. In this case, the graphene-adsorbate complexes always exhibit a 
spin moment of 1 /xg. The induced spin polarization texture is shown in Fig. [ljfor the 
adsorption of the groups adenine (b) and methyl (c), and the PMMA (d) and PTFE (e) 
polymers. The calculated patterns are remarkably similar in all the cases. Surprisingly, 
although the spin moment is induced by the functionalization, it is mostly localized in 
the graphene layer. The p z character of the electronic states behind the spin polarization 
is also evident in the plots. The distribution of the spin moment follows the bipartite 
character of the graphene lattice: C atoms in the opposite (same) sublattice than the 
saturated C atom Ci show majority (minority) spin polarization. Although the total 
spin moment is 1.0 fig, a Mulliken analysis only assigns ~0.10 /x# to the C atoms that 
form the graphene-adsorbate bond (Ci— X, with X=C). In contrast, the contribution 
from the three first nearest-neighbors [C2 in Fig. [TJ (a) ] is 0.34 /x#, -0.13 \1b from the 
next nearest-neighbors, 0.26 fiB from the third neighbors, and 0.40 \xb integrated over 
larger distances. This clearly shows the slow decay of the spin moment induced in 
graphene by this type of covalent functionalization. The latter is also reflected in the 
long-range exchange interactions between adsorbates that we describe below. 

In order to understand the origin of this spin polarization we have explored the band 
structure of the functionalized layers. Figure [TJ(f ) and (g) show the band structure for 
graphene functionalized with adenine and methyl groups, respectively, as representative 
cases. The magnetization comes from a very narrow defect state that is pinned at 
the Fermi level (E^). This state shows a predominant p z contribution from the nearest 
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neighbors C2 to the defect site, with a much smaller component from the adsorbate itself 
and a negligible contribution from the central carbon atom Ci. This state is reminiscent 
of the defect level that appears associated with a carbon vacancy in a 7r-tight-binding 
description of graphene [10]. The origin of the level is also similar in the present case, 
since the new C— C bond formed upon adsorption saturates the p z orbital at Ci. The 
small dispersion of this jt^-defect band combined with its partially filled character favors 
the spin- uncompensated solution. Tabled] shows the spin-splitting of the p z defect-band 
(5E S ) at T point [56]. This splitting varies in the small range 0.19-0.24 eV for all the 
adsorbates anchored by a single C— C bond (see also Fig. [2]). This further confirms the 
similar character and localization of the defect state behind the spin polarization for all 
the molecules. The energy gain with respect to the spin-compensated solution (AEj^ = 
Epm-Efm) is comparable for all adsorbates, being ~45 meV. 

The observed spin moment is localized in graphene and, therefore, it is derived 
from the electronic states of the carbon layer. However, the spin moment appears due 
to the covalent functionalization with non-magnetic molecules. As mentioned above, 
the adsorbate saturates the p^-orbital of one of the C atoms in the layer and, thus, 
creates a defect analogous to a 7r-vacancy. Still, according to the results in Table [TJ 
the polarity of the graphene-adsorbate bond is an additional ingredient that determines 
the size of the spin moment of the system |57J. As the ionic character of the bond 
increases (going from H to F in Table [1]) we observe a transition from a magnetic to a 
non-magnetic ground state. We can use the results from Mulliken population analyses 
to characterize the relative polarity of these bonds. Although the absolute values of 
the charge transfer have to be taken with care, since they are known to be strongly 
dependent on the size and quality of the basis set, results obtained with the same type 
of basis set can be compared and used to establish trends [58]. In our calculations 
both H and NH 2 radicals present almost negligible charge transfers of 0.003 and 0.046 
electrons from graphene, respectively. Consequently, H and NH2 behave similarly to 
the adsorbates presented previously, i.e., those bonded to graphene through homopolar 
C— C bonds. The electronic band-structures induced by H and NH2 are similar to 
those presented in Fig. [TJ and solutions with large spin-moments close to 1 fiB are also 
stabilized. In contrast, the charge transfers from graphene to OH and F increase to 0.26 
and 0.27 electrons, respectively. As a consequence, the Fermi level and the p 2 defect- 
band pinned to it appear now 0.2-0.3 eV below the Dirac point. This is consistent with 
the appreciable charge transfer towards the adsorbate and the consequent doping of 
the layer with holes. The additional population of the defect level, as well as its energy 
position in a region with a larger density of states of graphene, causes a drastic reduction 
of the spin moment for OH and a non-magnetic ground state for F. However, calculations 
using the fixed spin-moment method, presented in parentheses in Table [TJ show that 
the energy penalty needed to develop a spin solution is quite small even in the case of 
F (AEj\/=— 24.1 meV). Thus, although non-magnetic according to our calculation, the 
graphene— F bond can be easily polarized (at least in the low adsorbate-density regime 
explored here). This can be important to interpret the recent experimental report of 



Universal Magnetic Properties of sp 3 -type Defects in Covalently Functionalized Graphene! 



Table 1. Results of the structural parameters for all the studied adsorbates. We 
follow Fig. [T] (a) for the nomenclature. Energy gain AE^/ due to the spin polarization, 
total spin moment S and spin-splitting <5E S of the defect level are also included. The 
numbers in parentheses for F are calculated using the fixed-spin method. 



MA) S x (°) 8 2 (°) AE M (meV) S(/i B ) 5E s (eV) 

X = c 



CH 3 


1.59 


104.5 


113.9 


48.6 


1.00 


0.23 


C 2 H 5 


1.59 


105.0 


114.0 


48.0 


1.00 


0.23 


CeHn 


1.67 


105.9 


112.5 


47.4 


1.00 


0.23 


CeH 5 


1.59 


105.7 


112.9 


46.2 


1.00 


0.23 


CeH 4 F 


1.61 


106.6 


112.2 


45.6 


1.00 


0.23 


C 6 H 4 N0 2 


1.60 


107.0 


111.8 


32.3 


1.00 


0.20 


CgH40CH3 


1.59 


106.9 


111.9 


38.4 


1.00 


0.21 


C6H4CH3 


1.61 


106.6 


112.2 


46.6 


1.00 


0.23 


C 6 H 4 NH 2 


1.60 


106.5 


112.3 


42.8 


1.00 


0.22 


CONH 2 


1.65 


104.9 


113.6 


43.2 


1.00 


0.21 


COOH 


1.60 


104.5 


113.6 


41.6 


1.00 


0.22 


PMMA 


1.67 


106.8 


112.0 


47.1 


1.00 


0.23 


PTFE 


1.67 


105.8 


112.9 


64.8 


1.00 


0.19 


Adenine derivative* 


1.65 


104.7 


113.8 


42.5 


1.00 


0.23 


D-Glucose 


1.64 


105.1 


113.5 


41.4 


1.00 


0.21 


Polystyrene 


1.60 


106.8 


112.0 


40.0 


1.00 


0.21 


Polyacetylene 


1.58 


105.0 


113.6 


50.1 


1.00 


0.24 


X = H, N, 0, F 














H 


1.12 


102.6 


115.3 


46.5 


1.00 


0.24 


NH 2 


1.52 


105.2 


113.4 


27.7 


0.89 


0.20 


OH 


1.52 


103.7 


114.6 


8.4 


0.56 


0.12 


F 


1.55 


102.1 


115.7 


0.0 


0.00 


0.00 










(-24.1) 


(1.00) 


(0.16) 



* 9- (2-aminoethyl) adenine anchored to a carbonyl group [42J 



colossal magnetoresistance in this system. [59] 

We now look at the electronic structure around Ej? in more detail. Fig. [2] shows the 
calculated density of states (DOS) per spin channel for different adsorbates chemisorbed 
on graphene through a C— C bond. Despite the curves being shifted and smoothed with 
a Lorentzian broadening, the data collapse onto a single pattern: one fully polarized 
peak appears close to Ep- This confirms the universality of the spin moment associated 
with this type of covalent functionalization, independent of the particular adsorbate. 
In particular, our results point out to an almost perfect analogy between chemisorbed 
hydrogen [9j[T6l[T0] and adsorbates bind to graphene through single C— C bonds, which 
is not an obvious behaviour. Deviations (not shown in Fig. [2]) from this collapsed curve 
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Majority spin 




E-E F (eV) 

Figure 2. (Color online) Spin-polarized density of states for one molecule chemisorbed 
on a graphene supercell. All the adsorbates considered attach to the layer through a 
single C— C bond. For clarity, curves for different adsorbates have been shifted and 
smoothed with a Lorentzian broadening of 0.12 eV. The Fermi level is marked by the 
dashed line and is set to zero. The shaded regions denote the density of states of 
pristine graphene. 

are observed for adsorbates that form with graphene bonds of a considerable different 
character and polarity. 

Next we study the magnetic couplings between adsorbates at low concentrations. 
Figure [3] shows the energy differences between ferromagnetic (FM) and antiferromagnetic 
(AFM) alignments of the adsorbate-induced spin moments as a function of the distance 
between the adsorption sites. These results can be used to estimate the size and 
the distance dependence of the exchange couplings between the adsorbate-induced 
moments in graphene. We have studied systems of two different kinds: in one case 
two independent molecules are adsorbed at different locations, while in the other case 
a single molecule forms two covalent bonds with the graphene layer. To exemplify the 
first type of system, we have used H and CH 3 as representative adsorbates. Several 
observations can be made in this case: (i) If the two molecules are located at the same 
sublattice (AA adsorption), the FM alignment is always more stable than the AFM 
one. In the FM case, the total spin moment integrates to 2.00 [Xb for both H and CH3 
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Figure 3. (Color online) Schematic representation of some of the systems used to 
study the exchange coupling between adsorbate-induced spin moments in graphene: 
(a) two methyl groups, (b) the two ends of an alkane chain attached at two 
neighboring locations, and (c) the two bonds formed in a cycloaddition. The calculated 
energy difference, using a 8x8 supercell, between antiferromagnetic and ferromagnetic 
alignments of the local spin-moments are shown as a function of the distance along 
the armchair (d) and zigzag (e) directions. Only data for saturation at the the 
same graphene sublattice (AA configurations) are shown. Simultaneous adsorption 
at different sublattices (AB configurations) always gave rise to non-spin-polarized 
solutions. Filled and empty squares correspond, respectively, to H and CH3. Diamonds 
(green symbols) in panel (d) stand for the alkane chain and the star in panel (e) for 
the hypothetical AA-cycloaddition shown in (c). Circles (gray symbols) correspond to 
the best fit using a Heisenberg model (see text). 

and the local spin population remains nearly constant around each defect site, roughly 
insensitive to the presence of the neighboring adsorbate; (ii) If the two molecules are at 
the different sublattices (AB adsorption), we could not stabilize any magnetic solution, 
and the system converges to a spin-compensated solution with no local moment; (iii) 
The size and distance dependence of the magnetic couplings are very similar for H and 
CH 3 , which further supports the analogy presented above; (iv) For two methyl groups 
at short distances (< 5.0 A) we could only find FM solutions. 

We have also considered the magnetic couplings between local spin-moments 
induced in the layer by the formation of two single C— C bonds between graphene and 
the same molecule. We first explore this situation using a long alkane chain attached to 
graphene at two different locations [see Fig. El (b)]. Since the characteristics of the spin 
moment induced by covalent functionalization are quite independent of the adsorbed 
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molecules, as far as they are anchored by a single C-C bond, we expect in this case a 
result similar to that previously found for two methyl groups. This is confirmed by the 
results in Fig. [3] (d). For AA adsorption, i.e., when the two bonds between the alkane 
chain and graphene are established with atoms in the same graphene sublattice, the 
FM alignment with a 2 fig total spin moment is still the most stable. The size of the 
exchange couplings is almost identical to the case of two independent methyl groups. 
The only significant difference is that for the alkane chain we succeeded to find an AFM 
configuration at short distances. Furthermore, in the present case AB adsorption also 
converges to spin-compensated solutions. 

Figure [3] (c) and (e) show another example of magnetism associated with several 
bonds formed by the same molecule. In this case, we have considered a dipolar 
cycloaddition [60j|6l|. Usually this kind of cycloaddition saturates simultaneously two 
first-neighbor carbon atoms in graphene, i.e., it corresponds to an AB configuration. 
As a consequence, in agreement with the results presented so far, the cycloaddition 
does not create a spin moment in the layer. However, we have considered here another 
configuration in which the molecule is attached to two second-neighbor carbon atoms 
in graphene. This artificial AA-cycloaddition is 1.18 eV less stable than the standard 
AB-cycloaddition. In accordance to the universality of our results so far, the AA- 
cycloaddition gives rise to a spin-polarized solution. As expected, the FM solution with 
2fj,B spin moment is more stable, and the estimated exchange coupling fits well with 
those found for H, the methyl radical and the alkane chains. 

The observed magnetic behavior follows closely the expectations based on Lieb's 
theorem for a bipartite Hubbard model at half-filling [23], which is an appropriate model 
to describe the low energy electronic excitations in graphene. Lieb's theorem predicts 
that graphene should develop a total spin moment of 2/i# for AA adsorption and zero 
for AB adsorption [2U [TP] . However, this does not explain why a spin-compensated 
solution is more stable than an AFM solution for AB adsorption, as we find in our 
first-principles calculations. The reason for this behavior can be rationalized as follows. 
When adsorption takes place at opposite sublattices, due to the bipartite character 
of graphene, the interaction between defect levels in neighboring adsorption sites is 
appreciable. This interaction opens a bonding-antibonding gap in the p 2 -defect bands 
that contributes to stabilize the spin-compensated solution. The non-magnetic solution 
is favored when this gap is larger than the spin splitting 5E S of the isolated defect [62J. 
A detailed analysis of the calculated band structures confirms this interpretation. It is 
noteworthy that a similar behavior was also observed for chemisorbed hydrogen [9] [16] , 
vacancies [63j HQ], and even for substitutional Co atoms in graphene [62]. All of them 
exhibit a very similar electronic structure nearby E^. According to our interpretation, 
at sufficiently large distances between the adsorbates, a magnetic configuration with 
non-zero local spin-moments at each adsorption site should become favorable even for 
AB adsorption. At such long distances, the gap in the defect band becomes smaller 
than 5E S and the local spin-moments are recovered. However, in accordance with Lieb's 
theorem, an AFM configuration will be favored in this case. The necessary distances to 
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stabilize these AFM solutions are larger than those explicitly considered in our ab initio 
calculations presented in Fig. [3j The co-existence of non-magnetic and AFM solutions 
for AB configurations has been observed using first-principles calculations in the case of 
substitutional Co atoms in graphene [62J. 

An important consequence of the results presented above is that it is possible to 
turn graphene magnetic by chemical functionalization [21]. Unfortunately, they also 
tell us that, in order to achieve and control such magnetism, it is necessary to have 
selectivity on the adsorption site, i.e., molecules should preferentially attach to one 
of the graphene sublattices. This can be difficult to achieve in practice, particularly 
if we take into account that AA configurations are usually less favorable than AB 
configurations [161 EJ- Therefore, some difference must be introduced between A and 
B sites to promote adsorption on a single sublattice. For example, one possibility could 
be the registry with an appropriate substrate, but this can also alter other desired 
properties of the graphene layer. Fortunately, our results also show that it is not 
necessary to adsorb several molecules to obtain large spin moments in graphene. A 
similar effect is obtained by adsorbing a single molecule that simultaneously establishes 
several C— C covalent bonds with the layer. This simple observation might open a new 
route to engineer the magnetism of graphene derivatives. In principle, molecules could be 
synthesized with the appropriate distances between anchoring groups so, when adsorbed 
on graphene, all bonds will be formed with atoms in the same graphene sublattice. The 
spin moment associated with the adsorption of such molecules would be equal to the 
number of C— C bonds with graphene. Furthermore, the stability of this spin moment 
will be enhanced due to the strong FM interactions between nearby localized moments 
in the same sublattice. 

We have also analyzed our results for the magnetic interactions within the 
framework of a classical Heisenberg model. The interactions are described by the 
model Hamiltonian H = Yli<j <^( r jj)SjSj, where J{vij) is the exchange constant and 
Si is the local moment induced by covalent functionalization. The expression for the 
angular dependence of the exchange has been taken from the RKKY-like coupling 
obtained analytically in Ref. jH]. A simple |rjj| a behavior was assumed for the distance 
dependece, and the exponent a was fitted to our ab initio results in the case of AA- 
adsorption. We find that the interaction between the chemisorbed molecules in graphene 
surface is long range and falls off slowly with the distance, roughly proportional to 
JAA^ij) ~ l r ijl~ ( ' 1+e ' > with e ~ 0.20. [see the circles in Fig. [3] (d) and (e)]. This distance 
decay agrees well with a recent theoretical study using a 7r-tight-binding model [65J, 
where the interaction between adatoms on graphene was shown to decay as the inverse 
of the distance. In contrast, exchange interactions between substitutional Co impurities 
in graphene, a system that in principle can also be assimilated to a tt- vacancy or H 
adsorption, were recently shown to decay much faster (~ l^jl -2 ' 43 ) [62] • This is probably 
related to the larger localization of the spin moment in the latter case, associated with 
the 3d Co contribution. 
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4. Conclusions 

We have shown that covalent functionalization can be used to switch on the magnetism 
of graphene. Notwithstanding its structure and chemical or biological activity, when 
a molecule chemisorbs on graphene with the formation of a single C— C covalent bond 
between the adsorbate and the layer, a spin moment of 1.0 \ib is always induced in the 
system. This universal spin moment is almost exclusively localized in graphene, and its 
origin can be traced back to the polarization of a defect state that appears nearby the 
Fermi level. A similar effect accompanies the functionalization of carbon nanotubes [66J , 
and the chemisorption of hydrogen on graphene. This behavior is also reminiscent of 
the polarization observed for a single vacancy in a 7r-tight-binding description of the 
carbon monolayer. Other types of bonding, mediated by species different from C, also 
show similar magnetic properties as far as the established bond is not strongly polar. 
If chemisorption is accompanied by a substantial charge transfer between the adsorbate 
and graphene, the spin polarization is reduced and eventually disappears, as it is the case 
of F adsorption. However, even in this case, the low energy cost necessary to polarize 
the system seems to indicate that the magnetic properties of fluorinated graphene can 
be tuned using moderate magnetic fields. [59] 

Similar to the moment formation, exchange interactions are also found to be 
very weakly dependent on the nature of the adsorbates. In particular, we have 
explicitly checked that hydrogen atoms and methyl radicals give rise to the same 
magnetic couplings. Adsorbates in the same graphene sublattice (AA adsorption) 
couple ferromagnetically with an exchange interaction that decays almost inversely 
proportional to the distance JAA( r ij) ~ kiil^ 1 with e ~ 0.20. In contrast, the 
larger effective electronic hoppings between adsorbates at different sublattices (AB 
adsorption) prevent magnetism from arising in that situation. We have checked that 
the same results are obtained irrespective of whether the C— C bonds are established 
by the simultaneous adsorption of several molecules or by the adsorption of a single 
larger molecule. Therefore, the observed behavior is intrisic to graphene and appears 
associated with the saturation of a given collection of p z orbitals in the layer. This 
might open new routes to engineer the magnetism of graphene derivatives since, in 
principle, molecules can be produced with the appropriate structural and chemical 
characteristics to simultaneously bond to several atoms in the same graphene sublattice. 
The spin moment arising from such an adsorption will be proportional to the number of 
covalent bonds established and will have an enhanced stability due to the ferromagnetic 
interactions between nearby saturations in the same sublattice. Our findings are 
particularly attractive in the light of the recent sucessful synthesys of different graphene 
derivatives using surface chemical routes [30} 151) . Therefore, the synthesis of carbon 
nanostructures with functional groups at predefined positions, starting from previously 
functionalized monomers, seems plausible nowadays. 
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